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We show that preferential rewiring, which is supposed to mimick the behaviour of financial agents, 
changes a directed-network Ising ferromagnet with a single critical point into a model with robust 
critical behaviour. For the non-rewired random graph version, due to a constant number of out-links 
for each site, we write a simple mean-field-like equation describing the behaviour of magnetization; 
we argue that it is exact and support the claim with extensive Monte Carlo simulations. For the 
rewired version, this equation is obeyed only at low temperatures. At higher temperatures, rewiring 
leads to strong heterogeneities, which apparently invalidates mean-field arguments and induces 
large fluctuations and divergent susceptibility. Such behaviour is traced back to the formation of a 
relatively small core of agents which influence the entire system. 


I. INTRODUCTION 

To understand a complex behaviour of financial mar¬ 
kets is one of the main objectives of econophysics. Fat- 
tail non-Gaussian fluctuations, volatility clustering or 
rapid decay of autocorrelations of returns characterize 
most of the financial markets, suggesting that these styl¬ 
ized facts [l| have some more fundamental explanation. 
Searching for such an explanation, one can resort to the 
approach mrticularly suited for physicists, namely agent 
modeling [^. In the spirit of statistical mechanics, one 
considers a collection of agents involved in interactions 
resembling functioning of financial markets. Since buy¬ 
ing and selling are the most important activities of such 
agents, a number of models of financial markets bear 
some similarity to the two-state percolation Q or Ising- 
like models Q. 

An important agent’s interaction is mimicking some 
other agents’ behaviour, which suggests a similarity to 
ferromagnetic systems. However, ferromagnets typically 
exhibit rather small fluctuations, which is much differ¬ 
ent from the behaviour of finacial markets. Ferromag¬ 
nets exhibit large fluctuations only at the critical point 
separating ferromagnetic and paramagnetic phases. To 
place the system at the critical point requires, however, 
a fine tuning of control parameters. On the other hand, 
financial markets seem to be more robust with strong 
fluctuations appearing without any tuning of parame¬ 
ters. Some models were proposed, where mimicking other 
agents’ behaviour is compensated with the tendency to 
be in the minority Q or where agents with more complex 
strategies were used Q . They do reproduce some of the 
stylized facts but their considerable complexity hinders 
deeper understanding. 

Apparently, the analogy with simple ferromagnetic sys¬ 
tems is not sufficient to model financial markets and one 
should search for more suitable extensions. In our opin¬ 
ion, an important ingredient of models of financial mar¬ 
kets should be the possibility to choose and sometimes 
also change neighbours that a given agent would like to 
mimick. The objective of the present paper is to im¬ 
plement such a rewiring mechanism and to show that it 


drastically affects the behaviour of the model. When the 
neighbours to be mimicked are selected at random and 
kept fixed, the model behaves as an ordinary ferromagnet 
with ferromagnetic and paramagnetic phases separated 
at a critical point. However, when agents might switch 
the neighbours and preferentially select those they con¬ 
sider as more influential, the system generically exhibits 
divergent fluctuations. Such behaviour indicates that 
preferential rewiring induces a robust criticality, which 
is a required feature of stock-market models . We also 
examine the mechanism leading to the robust criticality. 

II. MODEL WITHOUT REWIRING 

In our model we consider N agents represented by spin¬ 
like variables Si = ±l, i = 1,2 ..., N. At each time step 
t, each agents decides whether to buy {si = 1) or sell 
(sj = — 1) an asset. To make the decision, an agent tries 
to mimick the behaviour of its neighbours and the model 
evolves according to the heat-bath dynamics: 

s ft + n = / ^ probability p = 

* 1 ~ 1 with probability 1 — p 

( 1 ) 

where 

hiit) = Y^Sj{t) ( 2 ) 
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is the local field acting on a given agent i and summation 
in Eq. m is over its neighbours. The control parameter T 
is the analogue of the temperature in the magnetic Ising 
model and determines the level of fluctuations in the de¬ 
cision process. 

The neighbourhood of a given agent is set randomly, 
namely, each agent has a fixed number of z randomly 
selected neighbours, which it interacts with via the lo¬ 
cal field. The neighbouring relation is not necessarily 
symmetric: if agent j enters the expression for the lo¬ 
cal field of agent i, it does not imply that agent i enters 
the expression for the local field of agent j. In other 
words, agents are nodes of a directed random network 
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and each node has z out-links (arrows point at the nodes 
that contribute to the local field). The number of in-links 
of a given agent, which specifies how many agents it in¬ 
fluences, is not fixed and it can vary among agents (of 
course, the average over all agents equals z). Equal num¬ 
bers of out-links and unequal numbers of in-links consti¬ 
tute an important feature of our model, which we will 
refer to as the out-homogeneity. 

Taking into account the spin variables, the above rules 
define actually an Ising ferromagnet on a directed ran¬ 
dom graph. Models of this kind were already analysed 
and shown to exibit an ordinary ferro-para phase transi¬ 
tion belonging to the mean-field universality class Q. 

In the following, we present a more detailed analysis 
of our model for z = A. Due to the out-homogeneity, one 
can write a relatively simple equation, which governs the 
evolution of magnetization. Let Pi{t) denote the prob¬ 
ability that agent i at time t takes the value = 1. 
Assuming that Pi(t) is spatially homogenous and does 
not depend on i, from the heat-bath rules we obtain that 

P(t+1) = V f- r-K - 

^ ^ l + exp[-4(/c-2)/r] 

( 3 ) 

Of course, Eq. can be easily rewritten in terms of 
magnetization {m{t) = 2P{t) — 1), which is common in 
Ising-model studies. In the steady-state limit (t —>• oo). 
Eq. ([3]) becomes a 4-th order polynomial equation, which 
can be easily solved numerically (and with some more 
effort even analytically). Moreover, the critical temper¬ 
ature Tc can be found using the standard procedure of 
expanding the t ^ oo limit of Eq. ([3]) in the vicinity of 
the critical point. Elementary calculations reveal that Tc 
obeys 


2 = tanh(4/rc)+ 2tanh(2/Tc). (4) 

The solution of the equation above can be written as 
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Tr = 


ln(i±^)’ 

' 1 — X > 


( 5 ) 


where 



1 - 5; 


II-k3v^ 


+ 



( 6 ) 


We thus obtain Tc ~ 3.08982, approximately. 

The factorized form of the probabilities suggests that 
Eq. ([3]) is nothing more than the mean-field equation for 
our model and thus it is only approximate. This would 
certainly be the case for undirected graphs, where neigh¬ 
bours j and k of agent i are strongly correlated (since i 
contributes to the local fields of both j and k). For undi¬ 
rected random graphs, some insight into the behaviour 
of the Ising model can obtained using a replica method 
0 or some recurrence relations based on the similarity of 
random graphs to Cayley trees [13. On the other hand, 
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FIG. 1. In a directed random graph, neighbours j and k of 
node i are not more correlated than any other two randomly 
selected nodes. 



FIG. 2. Temperature dependence of the magnetization m for 
the rewired and non-rewired models (2 = 4). The data are ob¬ 
tained from Monte Garlo simulations and are compared with 
the numerical solution of Eq. ®. Simulation and equilibra¬ 
tion times were equal to 10^ Monte Carlo steps. 


in directed networks even though j and k are neighbours 
of i, they are not more correlated than any other two ran¬ 
domly selected nodes (Fig.®. Since the graph is sparse, 
we expect that in the limit N ^ oo, such correlations are 
negligible, and consequently, the factorization in Eq. ® 
should be legitimate. 

Monte Carlo simulations of our model confirm the 
above analysis (Eig. Calculating the magnetization 
for N = 10^ and z = 4, we find it in a very good 
agreement with m = 2P(t = oo) — 1 obtained from the 
numerical solution of the steady-state limit of Eq. ®. 
For T = 3, we made much more extensive calculations 
(Fig.®. The linear extrapolation N ^ oo based on sim¬ 
ulations for < 3 • 10® gives m = 0.34723(2), which 
is in perfect agreement with m = 0.347225... obtained 
from the numerical solution of Eq. I®. In our opinion, 
such agreement strongly supports the claim that Eq. ® 
is exact (at least in the limit t —>■ oo). 


III. MODEL WITH REWIRING 

The model analysed in the previous section behaves 
similarly to some other Ising-like models with ferromag- 
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FIG. 3. Magnetization m as a function of the inverse of size 
calculated for the non-rewired model of size N = 10^, 3 • 10^, 
10®, 3-10®, 10®, and 3-10®, with z = 4 and T = 3. Simulation 
and equilibration times were equal to 10^ and 10“* Monte Carlo 
steps, respectively. In the limit N —>■ oo, a perfect agreement 
with the solution of Eq. can be seen. 


FIG. 4. Time dependence of magnetization m for the rewired 
and non-rewired models {z = 4, N — 10“^). Simulations were 
made for T = 4, which for the non-rewired model is deep in 
the paramagnetic phase. 


netic and paramagnetic phases separated at the critical 
point. Our primary motivation is to modify such ordi¬ 
nary ferromagnets so that they would resemble the be¬ 
haviour of financial markets, at least to some extent. We 
are particularly interested in supplanting a fine-tuned 
critical point with a more generic critical behaviour, 
which would exist in some, possibly large, temperature 
range. So far our agents make the decision to buy or sell 
based on the observation of their z neighbours, and the 
assignement of these neighbours is fixed during the entire 
evolution of the model. In the present section we mod¬ 
ify this rule and allow to change the neighbours. The 
rewiring we use is preferential: each agent has its sta¬ 
tus equal to the number of in-links that are (currently) 
attached to it. The selection of a new neighbour takes 
place with probability proportional to its status [HI. A 
single step of the dynamics of our model is thus defined 
as follows: 

• update spin variables Si {i = 1,2,, N) according 
to the heat-bath algorithm (HD- 

• rewire each agent selecting preferentially anew its 
z out-links. 

Since we keep the dynamics of spin variables basically 
unchanged, one might expect that Eq. (|31) still describes 
the behaviour of our model. Monte Carlo simulations 
show that to some extent this is indeed the case (Fig. [^l 
and a very good agreement with Eq. can be seen 
over much of the temperature range. However, close to 
the critical point T = T^, the rewired model shows much 
lower and perhaps zero magnetization. It would be desir¬ 
able to understand the reasons why Eq. m is no longer 
obeyed at higher temperatures. Possible explanations 
include appearance of correlations (that we argued are 


negligible in the non-rewired case) or a breakdown of ho¬ 
mogeneity (which is also one of the assumptions leading 
to Eq. ([3])). Some efforts to understand the origin of this 
behaviour will be made in the next section. 

What is even more interesting, the magnetization in 
the rewired version shows large fluctuations also at tem¬ 
peratures much higher than (Fig.lT]). To measure these 
fluctuations more quantitatively, we calculated the sus¬ 
ceptibility X that up to the temperature factor is equal 
to the variance of magnetization x = ~ 

5'i)^]- Numerical values indicate that as a function 
of system size N the susceptibility diverges as y ~ fV“, 
where a ^ 0.67 — 0.91 depends slightly on tempera¬ 
ture (Fig. [5]). Such behaviour is observed in a large 
temperature range (3 < T < 6) for the system size 
10^ < N < 3 - 10"^ . The divergence of susceptibility indi¬ 
cates that the model exhibits a robust critical behaviour. 
Together with data from Fig. [31 this suggests that the 
model with rewiring has two phases: low-temperature, 
which is ferromagnetic, and high-temperature, which is 
critical. It is difficult for us to locate precisely the transi¬ 
tion point between these two phases. For longer sim¬ 
ulations, it seems to shift slightly toward lower tem¬ 
peratures. Moreover, one cannot exclude that at suffi¬ 
ciently large temperature the critical phase will be re¬ 
placed with the paramagnetic one (having much smaller 
fluctuations). 

The critical behaviour in our model is also robust with 
respect to the frequency of rewiring. We made simula¬ 
tions with rewiring taking place, e.g., with probability 0.1 
(i.e., with probability 0.9, the out-links of a given agent 
at a given step were left unchanged). Such modification 
slows down the dynamics but retains the power-law di¬ 
vergence of the susceptibility. 
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FIG. 5. Size dependence of susceptibility x for the rewired 
model (2 = 4). The power-law £t x ~ N°‘ shows that a varies 
from 0.67 for T = 6 up to 0.91 for T = 3. 


IV. DYNAMICS OF REWIRING 

To get some understanding of our model, we looked 
at the structure of the network that emerges during the 
rewiring. Let us notice that rewiring is not affected by 
spin variables and thus might be considered as an inde¬ 
pendent process (but not vice versa—spin dynamics de¬ 
pends of course on the structure of the network). Some 
insight is already obtained from simulations of a small 
system (Fig. [5]). One can notice that most agents have no 
in-links and thus they do not influence any other agent. 
There is only a small core of agents which are responsi¬ 
ble for the decision formation of the other agents. Such 
structure appears also for larger systems (Fig. [T]). One 
can notice a substantial heterogeneity of the resulting 
core as for the number of agents that a given agent is 
influencing. 

A more detailed analysis shows, however, that the core 
size L slowly diminishes in time (Fig. [5]). This is not 
surprising since once an agent looses all of its in-links, 
it cannot get them back because the probability to be 
selected in a rewiring process is proportional to the cur¬ 
rent number of in-links (which is 0 for such an agent). 
Although the process of diminishing of L is irreversible, 
it is extremely slow for z > 1. Only at z = 1, this process 
is considerably faster and in a large time interval consis¬ 
tent with decay. Such a slow decay for z > 1 suggests 
that at long (but not infinitely long) time, the core is al¬ 
most in a steady state and has a certain size. Numerical 
calculations show that for z = 2 and 4 it increases with 
the system size approximately as (inset in Fig. 0). 

Some insight into the stability of the core can be ob¬ 
tained from the analysis of the time t needed for the 
system to condensate, i.e., for a given z to reach the core 
size z -I- 1 (which is the smallest core size that the system 
can reach). Numerical calculations show that for z > 1, 
T exhibits a fast, possibly exponential, increase with the 


FIG. 6. Network structure after simulations of t — 10® steps 
for A = 30 agents and z = 2 [l^ . Nodes with no in-links 
(open circles) represent agents which do not inluence decisions 
of any other agent (links that go out from them are drawn 
with dotted lines). Nodes with some in-links (filled red circles) 
represent the only agents that influence other agents (their 
out-links are drawn with solid lines). 



FIG. 7. Network structure after simulations of t = 10“* steps 
for N = 10® agents and z = 2 [13] . Agents with no in-links 
are omitted. The size of a circle is proportional to the number 
of in-links. 


system size (Fig. [9]). And again, a slower increase (^ N) 
is obtained only at z = 1. 

Our results in Figs. [HHSlshow that rewiring for z = 1 
leads to a rather fast condensation while for z > 1 the 
dynamics is basically trapped in a core of size ^ 

Even though the condensed state could be reached in 
principle, for large JV and z > 1 it virtually never hap¬ 
pens. The situation is reminiscent of some models with 
the so-called absorbing states: for some values of con¬ 
trol parameters, the absorbing state of the dynamics is 
basically unavailable and the model remains in the ac¬ 
tive phase (the lifetime of which in that regime is also 
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FIG. 8. Time dependence of the size of the core (log-log 
scale); calculations were made for N = 10"^. Inset shows that 
for 2 = 2 and 4 the ’steady-state’ size of the core increases 
approximately as To generate the network of fractional 

2 (1 < 2 < 2), with probability 2 — 2 we created one out-link 
and with probability 2 — 1 we created two such links for each 
agent. Thus fractional 2 has only an average sense. 



logio(N) 

FIG. 9. The system size dependence of the time r to conden¬ 
sate, ie., to reach the core size 2-1-1. Let us notice that even 
for 2 slightly larger than 1, r shows a rapid increase. 

exponentially divergent with the system size) [l3j |. 

A network structure with a nearly stable core suggests 
an explanation of the generic divergence of susceptibility 
that we reported in the previous section. Indeed, since 
agents are influenced only by agents from the core, on av¬ 
erage there are agents that are influenced by a sin¬ 
gle agent belonging to the core. Considering core agents 
as independent (and influencing other agents), we 

easily obtain that y ^ Numerical results (Fig. [5]) 

suggest a faster increase at low temperatures (for T = 3, 
we obtained x ^ that most likely result from some 

correlations between core agents. Another factor affect¬ 
ing our simple estimations of the divergence of y might 
be some heterogeneity of the core (Fig. [7]). 


V. CONCLUSIONS AND REMARKS 

In the present paper we have shown that preferential 
rewiring changes an Ising ferromagnet, which has a sin¬ 
gle critical point, into a model with robust critical be¬ 
haviour. The rewiring mechanism that we used is sup¬ 
posed to mimick the behaviour of financial agents who try 
to follow their neighbours but at the same time have also 
some freedom to choose the ones to follow. We assume 
that the preference in the rewiring process is proportional 
to the number of in-links of a given agent. It is thus not 
a (more or less) objective measure of its performance but 
solely how the agent is perceived by the population of 
other agents. Similar recipes turned out to be successful 
im e.g., some page rank algoritms used by search engines 
[T3 | or various recommendation systems [l^. 

Our model is of some interest from the statistical- 
mechanics point of view. Due to out-homogeneity, we 
could write a simple mean-field-like equation (Eq. (|5|) ) 
that can be used to obtain the magnetization of the 
model. We argued, however, that for the present model 
in the non-rewired version this equation should be exact, 
and numerical simulations provide a very strong support 
for the claim. A very good agreement with this equa¬ 
tion was obtained also for the rewired case, but only in 
the low-temperature regime. We made some attempts to 
explain why rewiring invalidates this equation at higher 
temperatures and at the same time leads to the diver¬ 
gence of suceptibility and criticality. In our opinion this 
is related with the formation of a relatively small sub¬ 
set of agents that retain some in-links and thus drive the 
entire system. 

The change of the dynamical regime in the rewired 
version at 2 = 1 is also of some interest. It is tempting 
to associate the change with some percolation transition 
that for random graphs is known to take place at 2 = 1 
[T^ . However, rewiring redistributes links in a highly 
nonrandom fashion and a possible relation with random 
graphs is by no means obvious. 

We hope that our model might be useful also in the 
econophysics context. Relatively simple rules that gen¬ 
erate a robust criticality might serve as a starting point 
for further modifications and analysis. For example, one 
might consider a model where an agent that no longer 
has any in-links still retains some (small) status e and 
can be thus selected during the rewiring process. We 
expect that for small (possibly A-dependent) e such a 
model would be similar to our (e = 0) model, but cer¬ 
tainly numerical simulations would be needed to support 
such a claim. One of the important stylized facts that ap¬ 
parently is missing in our model is volatility clustering. 
One might hope that some extensions where agents, for 
example, try to be in the minority (like in the so-called 
minority games [13) or use more sophisticated strategies 
(like ’fundamentalist’, ’trend follower’ or ’noise trader’) 
will provide a more realistic description of financial mar¬ 
kets and at the same time will retain simplicity of the 
model . 
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